Formal exact operator solutions to nonlinear differential equations 



YU. N. KOSOVTSOV 

Lviv Radio Engineering Research Institute, 79060 Lviv, Ukraine 
email: yunkosovtsov@gniail.com 



Abstract 

The compact explicit expressions for formal exact operator solutions 
to Cauchy problem for sufficiently general systems of nonlinear differential 
equations (ODEs and PDEs) in the form of chronological operator expo- 
nents are given. The variant of exact solutions in the form of ordinary 
(without chronologization) operator exponents are proposed. 



1 Basis of the method. 

The starting point for solving differential equations by the operator method is 
linear first-order differential equation for operator E = E(t): 

■^-um (1) 

with the initial condition 

nt^a = I, (2) 

where L(<) is a linear operator which, generally speaking, depends on t, but 
does not depend on operator ^ explicitly; I is identity operator. The solution 
of this equation for t > a was obtained in 1949 by F. Dyson in the chronological 
exponential form Ij (see also [2], [3]) 



E(t) = Texp{/" dTL(T)}, 

J a 



(3) 



where it is supposed, that exponent represents a power series expansion with 
following chronologization according to the rule 

T{L(ti) L(T2)...L(T„)}=L(r„J L(t„J . . . L(t„J. 

Very important role plays the equation adjoint to fl]) for operator = 
E-i(i) (see H, 0) 

2^ - -E-LW (4) 
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with the initial condition 

E-'\t^a=l- (5) 

It is easy to see using the Dyson's idea, that its solution is 



E-^(t) 



Toexp{- / drL(r)}, (6) 

J a 



where operator Tq is the operator of opposed chronologization with action 
To{L(ti) L(t2)...L(t„)} = L(t„J L(t„J...L(t,J. 

It appears that the operator E~^(t) is inverse to the E(t), i.e. 

Toexp{- / drL(T)}Texp{ / dTL{T)}=I 

J a J a 

Ordering operators T and Tq for t < a change places. In some sense in linear 
problems operators T and Tq provide description of direct and back scattering 
accordingly. It be assumed henceforth that t > a. 

Ordering (chronological) operators have to be introduced because the oper- 
ators in the integrands with different values of the variable r may not commute, 
what cause certain difficulties in manipulation with chronological exponents. 
Nevertheless, examination of algebraical properties of the chronological expo- 
nents allows to obtain an extensive family of operator identities |4j-[B] which 
provide possibilities of transformations and analysis of chronological expres- 
sions. 

The operator A is derivative if 

A(A + B)/ = AA/ + AB/, 

AA/ = (AA)/ + AA/ 

for any differentiable function / and any differentiable linear operators A and 
B. In other words, the linear operation is derivative if it satisfies to Leibnitz 
rule. As we will see later, the chronological operators with derivative operators 
in exponent, which we will denote for short as 

E = Texp{/" drA(T)}, 

J a 

play the extremely important role here. 

Let b{t) be a, function that is the operator of multiplication on the b(t). Then, 
as can be seen from the definition of the derivative operator, the commutator 
[b{t), A(r)] is a function too. 

Now consider the following construction 

K = Eb{t)'E~\ 



2 



As far as the commutator [b{t), A(t)] is a function, then all repeated commuta- 
tors are functions too. So with help of BCH formula (see (27) in [6]) we conclude 
that operator K is a function and 

K = (E6(i)), (7) 

where the outward brackets in the right-hand side stress that action of operator 
is restricted by the brackets area. An analogous conclusion holds for E^^ b{t) E. 
Let us consider further the following obvious chain 

E &i ^2 ■ • • fen = E 6i E"^ E 62 E"^ . . . E b„ = 
= (E6i)(E62) ... (E&„). 

With what is observed here and the fact that the chronological operator is 
linear, leads to conclusion that for any function F(6i, 621 ■ ■ • 1 &n)j which can 
be expanded in power series with respect to &i, h^^ . . . , 6„, we can obtain the 
following nice property: 

E^^(&i,62, 6„) =F((Efoi), (E62), (E6„)). (8) 
Analogously 

E-i i^(6i, 62, ... , 5„) = F{{^-^ 61), (E-i 62), . . . , (E-i 6„))- (9) 

This property of the operator E coincides with familiar property of the 
shift operator exp{s^}, as shift operator is the particular case of chronological 
exponent. 

It be assumed henceforth that Fi are analytic functions on all their argu- 
ments in corresponding regions. 

2 Chronological solutions of nonlinear differen- 
tial equations. 

In existent literature, devoted to operator approach for solving differential equa- 
tions, as a rule, the linear problems are considered, i.e., one proceed from ([1])- 
([3]). The connection of equations ([T]) and ([3]) with nonlinear ordinary differential 
equations seemingly for the first time was explicitly mentioned in [4], [5]. The 
similar attempt to nonlinear partial differential equations in [7] had not led to 
compact form of solution, only infinite chronological iteration series was ob- 
tained there. 

Let us consider particular cases of the derivative operators A and differential 
equations, which correspond them. 

2.1 Ai(r) = F(r,c)|:. 
Differentiating function 

u(<,c) = Toexp{- /" dTF{T,c)^}c (10) 
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on t, we obtain, taking into account ^ and then that 

— = To exp{- J dT F{t, c)-g^}F{t, c) = u), 



i.e, the expression (jlOp is a solution of ordinary differential equation of the first 
order for u = c) (see [1], [S], and also 

— + F(i,^.) = 0, (11) 



where c is an arbitrary constant, u\t=a = c. That is, the expression (jlOp is exact 
general operator solution of the equation (fTTj) . 

2.2 A2(r) = E;=ii^.(r,ci,...,c.)^. 

Analogously, differentiating functions u,; = Ui(i, ci, . . . , c„), where i = 1, . . . , n, 
Mj = Toexp{-/ dr ^Fj(T, ci,...,c„) — }ci (12) 



on t, we ascertain that (see 6 ) is the solution of the following system of 
ordinary differential equations of the first order for ui = Ui(t, ci, . . . , c„) 

^+F,(i,ui,...,M„) = 0, (i = l,...,n) (13) 
where Ci are arbitrary constants, Ui\t=a — 

2.3 As{r) = d-C E;=i FAt, C, c, . . . , D,"c,, . . . )4^, 

where t is dedicated independent variable (" time" variable) , — Ci, ■ ■ ■ , Cm are 
"space" variables, the functions c = c(C) = ci(C), ■ • • , c„(C) will be later on play 
role of arbitrary functions of "space" variables (it is clear, that it is supposed 
that all these functions are arbitrary smooth and have not singularities), the 
symbol . . . , D"ck, ■ ■ ■ denotes a set of all partial derivatives of {k = 1, . . . ,n) 
on " space" variables Q of the type 



IS tl 

Differentiating on t functions 



up to some certain order, — is the functional derivative. 

ocj (C) 



Ui{t, x) = To exp{- / dr A3} Ci{x) , 

^ a 



(14) 



4 



where x = xi, . . . , x^, we obtain, taking into account (U) and then ((S]), that 
-^ = -Toexp{-y drAajo 

/OO 
<rC 6{C~ x) F.{t, C, c, . . . , D^ck, ...) = 
-OO 

= -Toexp{- f dT A3} F4t,x,c,...,D2ck,... ) = 

J a 

= -Fi{t,x,u,...,D'l^Uk,...) 

and ascertain that (|14p is the solution of system of partial differential equations 
of the first order on variable t (and arbitrary order on "space" variables) for 

Ui(t,x) 

011 ■ 

-^+F,{t,x,u,...,D^Uk,...)^0, {i = l,...,n) (15) 

where u = ui{t, x), . . . , u„{t, x) and now in . . . , D'^Uk, ■ ■ ■ appears the corre- 
sponding set of partial derivatives of Uk{t, x) on "space" variables x. 

For equations of higher orders on t by substitutions of type = Vi very 
often one can get over to system of the first order equations and, if it possible to 
resolve such system with respect to derivatives of unknown functions (including 
Vi) on t (i.e., reducing the equation or system into the so-called normal form), 
it is easy enough by using 2.2 or 2.3 to obtain operator solutions for these cases 
too. 

The existence of formal solutions of nonlinear differential equations in the 
form of chronological operator exponents is important as a practical matter, 
first of all, because they represent in compact and clear form the dependence of 
solutions from all parameters of the problem. If given differential equation fits 
the conditions on existence and uniqueness of solutions, then all forms of exact 
solutions are equivalent (with minor reservations) with each other. However, in 
usage of one or another form of solution one is wishing to have the opportunity 
to analyse it and certain freedom in its transformation in follow-on calculations. 
The interpretation of chronological solutions as infinite series is the first, but 
far from being an only opportunity. There are simple examples, when operator 
solution can be transformed to the conventional form just by operator methods 
(see, e.g., [B]). 

Fortunately for the stage of analytic (symbolic) calculations there is an ex- 
tensive family of operator identities [l]-[6], which allow to carry out not only 
basic algebraic operations with chronological exponents, but accomplish its dif- 
ferentiation on parameters and so on. If you wish, algebraic transformations here 
can be interpreted as resumniation of corresponding infinite iteration series in 
which compactness, visibility are kept and, what is important, opportunity to 
remain as soon as possible in terms of exact solutions. 

Nevertheless on the stage of the decision- result ("number") obtaining very 
often it is necessary to turn to series. The presence of chronologization, which is 
overburdened by multiple integration, complicates considerably as calculations 
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as essential analysis of series convergence conditions. Therefore during more 
than half a century the attempts to represent operator solutions in form of 
ordinary (without chronological operator) exponent are made (see survey [5]), 
here one of typical approaches is the so-called Magnus expansion. There are 
some variants to solve this problem, which have its own benefits depending on 
posed a goal. As a rule such approaches lead to approximate presentation in 
the form of very bulky expressions. The variant proposed below is differed by 
noticeable simpleness and is, what is more, the exact solution. 



3 The solutions of nonlinear differential equa- 
tions in form of ordinary (without chronolog- 
ical operator) exponent. 



It is follows from operator identity (23) of [5] and well-known properties of shift 

d_ 
dt 



operator that (note, operator L(t) does not depend on operator 4r explicitly) 



from here 



exp{{t-a) — }exp{{t-a)[L{s) - —]} 



= T exp{ f dr L(r + s - a)} , 

J a 



exp{(t-a)[L(.)-^]} 



d /■* 

:exp{-(t - a) — }Texp{ / drLiT + s-a)} 

OS Ja 



/■* d 
=Texp{/ dTL{T + s-t)}exp{-{t~a) — }. 

J a OS 

or equivalent 

d d 
exp{(i-a)[L(s) - — ]} exp{(i - a) — } = 
OS OS 

= Texp{f drL(r + s-t)}. 

J a 

Similarly for inverse operators 
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d d 
exp{-(i-a)[L(s) - — ]} exp{-(t - a) — } 



Toexp{- f drL 

J a 



(r + s-a)}. (16) 

Let us apply the obtained identities to the solutions of differential equations 
of Section 2. Note, first of all, that if both sides of identity (fTB]) multiply from 
the right on constant Ck or function Ck{x), which do not depend on s, then the 
action of right shift operator in the left-hand side will be equal to the action of 
identity operator /, while left-hand side operator under substitution s = a turns 
into above examined operator ([S]). 

Thus for 

L(r) = Ai(r)=i^(r,c)^ 
we have, taking into account (|16p and foregoing remarks, that 



/■* d 
u{t) = [Toexp{- / dr F{t + s - a,c) — }c\\, 

J a 



d 

d d 

= exp{-(t - a)[F{s,c)— - — ]}c|s=a- 
oc OS 

As a result we get the solution of equation (fTTI) in exponential form (without 
chronological operator) 

uit) =exp{-(i-a)[i='(s,c)— - —]} c\s=a- 

oc OS 

By the same way one can obtain the exponential form of solution for system 
of ordinary differential equations (|13p 

" d d 

Ui{t) = exp{-(t - a)[^i^j(T, ci, . . • 'C„)^ - |s=a, (i = 1, . . . ,n). 

j=i ■' 

Similarly the exponential form of solution for system of partial differential 
equations ([TS]) is as follows 

d 

Ui(i, f) = exp{-(t - a)[A3(s) - — ]}Q(f) \s^a, (i = 1, . . . , ti). 

If expand the operator exponent into series, then one can obtain the expan- 
sion of the solution into Taylor series on powers of [t — a). This is most obvious 
(but not unique) interpretation of the operator exponential form of solutions. 

The Maple procedures for analytical calculations of approximate solutions of 
ordinary and partial differential equations (Cauchy problem), which are based 
on the approach considered here, are presented in 
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